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Abstract 



Ordinary orthogonal polynomials are uniquely characterized by the three term re- 
currence relations up to an overall multiplicative constant. We show that the newly 
discovered M-indexed orthogonal polynomials satisfy 3 + 2M term recurrence relations 
with non-trivial initial data of the lowest M + 1 members. These include the multi- 
indexed orthogonal polynomials of Laguerre, Jacobi, Wilson and Askey- Wilson types, 
^vq ■ The M = case is the corresponding classical orthogonal polynomials. 
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Exactly solvable quantum mechanical systems in one dimension had been well studied [1] 
[3]. In recent years the subject saw remarkable developments [I] [H|. The eigenf unctions 
of the solvably deformed systems are described by new types of orthogonal polynomials, 

d • 

the exceptional and multi-indexed orthogonal polynomials. The first breakthrough was the 
discovery of the X\ Laguerre and Jacobi exceptional orthogonal polynomials in the context 
of Strum-Liouville theory by Gomez-Ullate, Kamran and Milson jl]. Quesne constructed 
shape- invariant quantum mechanical systems whose eigenfunctions are described by the X\ 
Laguerre and Jacobi polynomials [5]. The second progress was the construction of the Xp 
Laguerre and Jacobi exceptional orthogonal polynomials for any positive integer £ by Sasaki 
and the present author [6] , which were based on the deformations of the quantum mechanical 
systems preserving the shape- invariance [2]. The third development was the generalization 
of the exceptional orthogonal polynomials, i.e., multi-indexed orthogonal polynomials of 
Laguerre and Jacobi types, which were obtained based on the method of virtual states 



deletion for quantum mechanical systems |8]. The exceptional orthogonal polynomial (in 
the narrow sense) is a one-indexed orthogonal polynomial. Parallel to the ordinary quantum 
mechanical systems, the discrete quantum mechanical systems had been developed [34] [36] 
and the exceptional and multi-indexed orthogonal polynomials of Wilson, Askey- Wilson, 
Racah and Racah types were constructed [3Tj [4~T]. 

By the method of virtual states deletion [8], which is based on the Darboux-Crum trans- 
formation [4~2| |4~3] . infinitely many exactly solvable quantum mechanical systems are system- 
atically obtained from the original exactly solvable systems. The multi-indexed orthogonal 
polynomials of Laguerre, Jacobi, Wilson, Askey- Wilson, Racah and g-Racah types describe 
their eigenfunctions. They satisfy second order differential or difference equations and form 
complete basis but their degrees start at a certain positive integer £ instead of zero, namely 
the set of degrees is {£, £ + 1,£ + 2, . . .} (a maximum degree exists for ((/-)Racah cases). 
Thus the constraints of Bochner's theorem [H] are avoided. The Krein-Adler transforma- 
tion (eigenstates deletion based on the Darboux-Crum transformation) also gives infinitely 
many exactly solvable quantum mechanical systems from the original exactly solvable sys- 
tems [1S1 [35] . If the original system is described by (ordinary) orthogonal polynomials, the 
eigenfunctions of the deformed systems give new orthogonal polynomials. They also satisfy 
second order differential or difference equations and form complete basis. The main differ- 
ence from the multi-indexed orthogonal polynomials is their degrees. The set of degrees of 
these new polynomials is {0, 1,2,.. .}\{a?i, d 2 , ■ ■ ■ , (1m}, namely there are several 'holes'. The 
shape-invariance of the original systems is lost when deformed by the Krein-Adler trans- 
formations, whereas the system retains shape-invariance when deformed by the method of 
virtual states deletion. Since the new orthogonal polynomials obtained by the Krein-Adler 
transformation are not so natural in these senses, we are mainly interested in the multi- 
indexed orthogonal polynomials. 

Some properties of the exceptional and multi-indexed orthogonal polynomials have been 
studied [U EH EJ E], [HI] P2]> However, there remain various properties to 

be clarified. In this paper we focus on the recurrence relations. The ordinary orthogonal 
polynomials ('0-indexed' orthogonal polynomials) are completely characterized by the three 
term recurrence relations [46] . The three term recurrence relations are basic properties of the 
orthogonal polynomials and it is important to find the corresponding recurrence relations 
of the multi-indexed orthogonal polynomials. We will show that the M-indexed orthogonal 
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polynomials of Laguerre, Jacobi, Wilson and Askey- Wilson types satisfy the 3 + 2M term 
recurrence relations (j3.4|) . The coefficients of the 3 + 2M term recurrence relations are de- 
termined in terms of those of the three term recurrence relations of the original polynomials. 
In order to obtain the whole multi-indexed orthogonal polynomials by using these 3 + 2M 
recurrence relations, we have to specify the first M+l members of the polynomials as inputs, 
which are severely constrained. It is expected that similar recurrence relations hold for the 
multi- indexed Racah and g-Racah polynomials |40j . 

This paper is organized as follows. The essence of the Darboux-Crum approach to quan- 
tum mechanical systems and the multi-indexed orthogonal polynomials is recapitulated in 
section [2j Section [3] is the main part of the paper. After recalling the three term recurrence 
relations for the Laguerre, Jacobi, Wilson and Askey- Wilson orthogonal polynomials, we 
present the 3 + 2M term recurrence relations of the multi-indexed orthogonal polynomials. 
The initial data for the recurrence relations are discussed in § 13.41 The final section is for a 
summary and comments. The explicit formulas of the multi-indexed orthogonal polynomials 
of Laguerre, Jacobi, Wilson and Askey- Wilson types are presented in Appendix. 

2 Quantum Mechanical Systems and Multi-Indexed 
Orthogonal Polynomials 

Not only the construction of the multi-indexed orthogonal polynomials but also the derivation 
of the recurrence relations are based on the quantum mechanical formulation. See, for 
example, [36] for the general introduction. Here we recapitulate the Darboux-Crum approach 
to quantum mechanical systems with a continuous dynamical variable x and the multi- 
indexed orthogonal polynomials |I2| |j3| l4"5| 134] l35| l8| l4T] . 

We consider quantum mechanical systems in one dimension (x\ < x < X2), 

H = A j A, H(j) n {x) =£nMx) (" e Z> ), = £ < £ l < S 2 < ■ ■ ■ , (2.1) 

d f f X2 

((j>n, <t>m) = / dx(f) n (x)*(f) m (x) = h n 5 nm (h n >0). (2.2) 

Any solution of the Schrodinger equation T-L(j)(x) = S(j)(x), which need not be square inte- 
grable, can be used as a seed solution for a Darboux-Crum transformation. Let us take M 
distinct seed solutions {(p v (x)}, 

H(j) v (x) = £ v <f) v (x) (v = dx, d 2 , . . . , d M )- (2.3) 
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The s-step Darboux-Crum transformation with seed solutions V (v = di, (1%, . . . , d s ) gives 
^di...d s = f Ad- L ...d s $d x ...d s +£d a , (2.4) 

<j>di...d a n( x ) — •Ad 1 ...d s 4>d 1 ...d s _ 1 n( x )> <f ) d 1 ...d s v( x ) = •Ad 1 ...d a 4'd 1 ...d a -i v(%) j (2-5) 
n d 1 ...d 3 4 > d\...d a n{.%) — £n < Pdi...d s n 

(x), T-Ld 1 ...d a ( l>di...d s v{%) — £v<i>di...d a v(x)- (2-6) 

Here the concrete forms of A dl ,,, ds are given in (12.141) and (12.211) . and the eigenfunctions 
<f>di...d s n(%) an d the seed solutions <f>di...d a v(%) are expressed by using the Wronskians ( 12.151) 
or Casoratians (I2.23H . Since these (12.4D — (12.6p are shown in algebraic way, ( 12.61) holds for any 
range of the coupling constants contained in the system. However, the Hamiltonian 1-Ld x ...d a 
may be singular in general. By picking up another seed solution <p ds+1 , the Hamiltonian 
T-tdu-ds is rewritten as H dl ... ds = A ] di ,„ da+1 Adi...d a+l +S ds+1 . After M steps, we obtain (H [M] = 
W-dL.-dM, 4> [ n\x) = <f) dl ... dM n(x)), 

H^\x)=£ n ^\x). (2.7) 
If this Hamiltonian is non-singular, we have 

M 

(0f ], 4>M) = ]J(£ n -£ dj )- h n 5 nm , (2.8) 

and the Hamiltonian can be rewritten in the standard form 7i dl ... dM = A dl dM A dl ... dM - 
Note that the deformed systems are independent of the orders of deletions (4> dl ... ds n(x) and 
0<ii...d s v(^) may change sign). 

We know several methods for constructing non-singular Hamiltonian W- M '. Here we 
mention two methods, (i) [Krein-Adler transformation] [4"5| [35] . As seed solutions, the 
eigenfunctions are chosen, 4> Y (x) = <fi v (x), £ Y = £ v , and the index set {di, . . . , c?m} is required 
to satisfy the Krein-Adler conditions rij=i( m — dj) > (Vm G Z>o). In this case the 
intermediate Hamiltonians H dl ...d a may be singular but the final Hamiltonian w- 1 ^ is non- 
singular. The eigenfunctions are <p\^\x) with n G Z> \{(ii, . . . , c/m}- Compared to the 
original system H, M states with energy £^ are missing in the deformed system 1-fi M \ 
(ii) [method of virtual states deletion] [El HI]. As seed solutions, the virtual state wavefunc- 
tions are taken. For the definition of the virtual states, see [8j |JT]. The Hamiltonian "H' M l is 
non-singular (The parameter range may be restricted.) and the eigenfunctions are </>L M '(a;) 
with n G Z> . The deformed system W- 1 ^ is exactly iso-spectral to the original system %. 
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In this case the intermediate Hamiltonians / Hd 1 ...d s are also non-singular and iso-spectral to 
the original system. 

Let us assume that the eigenfunctions of the original systems in § 12.11 and § 12.21 are of 
polynomial type: 

<t> n {x) = Mx)Pn{v(x)) {neZ>o)- (2.9) 

Here P n {v) is a polynomial of degree n in 77 and 77 = rj(x) is a certain function of x, which 
is called the sinusoidal coordinate [47J. Then the eigenfunctions of the deformed system are 
also of polynomial type, 

<l>M(x) = *M(x)pM( V (x)), (2.10) 
where P^\r]) == Pd^.^M^iv) is a polynomial in 77. They are orthogonal to each other, 

M 

dx^ M \x) 2 P^{r,(x))P^{r,(x)) = H(S n -£ dj ) ■ hj nm . (2.11) 

The degree of Pn^(rj) is different from n. For (i), the degree of Pn M \i]) is generically 
£ + n, where £ = YljLi dj ~ \M(M + 1). The label n does not take all values in Z>o- It 
takes n 6 Z> \{di, . . . , (1m}, namely there are M 'holes'. For (ii), the degree of Pl M \r]) is 
generically £ + n, where a positive integer £ is determined by {di, . . . , c^u}, ( 1A.4I) . The label 
n takes all values in Z>o and there are no 'holes'. For the original systems described by the 
Laguerre, Jacobi, Wilson and Askey- Wilson polynomials, we call the obtained polynomials 
the multi-indexed orthogonal polynomials for Laguerre, Jacobi, Wilson and Askey- Wilson 
types. The eigenfunctions ^>di...d a n( x ) °f the intermediate Hamiltonians 'Hd 1 ...d a also have the 
following form: 

(pd^dsn^x) = ^ dl ... da (x)P dl ...d a ,n{ri(x)) (nGZ> ). (2.12) 

In the rest of the paper we consider the method (ii) only. The quantum systems to be con- 
sidered have some parameters (coupling constants), denoted symbolically by A = (Aj., A 2 , . . .), 
4> n (x) = (f) n (x;X), P n (rj) = P n (r];X), etc. As a quantum mechanical system, the parameter 
range should be chosen such that the deformed system is non-singular. In this pa- 

per, however, we treat the algebraic aspects of the multi-indexed orthogonal polynomials 
and various algebraic relations hold independently of the parameter ranges. So we do not 
care much about the parameter ranges. The explicit forms of the multi-indexed orthogonal 
polynomials of Laguerre, Jacobi, Wilson and Askey- Wilson types are relegated to Appendix, 
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since they are somewhat lengthy and they are not directly needed for the derivation of the 
recurrence relations. 



2.1 Ordinary quantum mechanics 

For ordinary quantum mechanics, the Hamiltonian with zero ground state energy can be 
expressed in a factorized form: 

= AU, V=^-^log|0 o (x)|, = --^-d x log\Mx)\- (2-13) 
The intertwining operators A^..^ and A dl ds are given by 

A*!...*, = - d x log\4> dl ... ds (x)\, A\ i dg = - ^logl^...^^)!, (2.14) 

and the eigenfunctions <pd 1 ...d s n(x) and the seed solutions </>d 1 ...d 3 v(x) have fractional expres- 
sions 

, / \ W[0 dl ,...,0 ds ,0„](x) ~ W[0 dl ,...,0 da ,0 v ](x) 
<l>du..d.n{x) = — — = , </> dl ...d.v(x) = — — j ■ (2.15) 

W[0 dl , . . . , d J (x) W[0 dl , . . . , d J (x) 

Here W[/i, / 2 , . . . , f n ](x) is the Wronskian 

W[/ 1 ,/ 2 ,...,/ n ](x) d ^ f det(^f^) , (2.16) 

V dxl 1 / l<j,k<n 

and W[-](x) = 1 for n — 0. 

We consider two shape-invariant systems, the radial oscillator and Darboux-Poschl- Teller 
potential, whose eigenfunctions are described by the Laguerre (L) and Jacobi (J) polynomi- 
als. Various data of these systems are: 

L: < x < oo, A = g, 8 = 1, <? > -, 

U(x;X) = x 2 + ^^--(l + 2g), S n {\) = An, n{x) = x 2 , 
x A 

0o(x; A) = e"^V, P n ( V ; A) = L ( r h \v), 

r(n + </+§), (2.17) 

7T 1 

J: 0<x<-, \=(g,h), 8 = (1,1), g,h>-, 

U{x; A) = 9 -4 7 ~ + — 7 ~ ~ (9 + ^ £„(A) =4n(n + # + /i), ??(x) = cos2x, 
sin x cos z x 



M*; = (sinx)»(cosi)*, P„(ij; A) = rf s *■* *>(ij), 

= r ( , i + g + i)r(, i + ft + i) 

" v ; 2n!(2n + ^ + /i)r(n + ^ + /i)' v ; 

where Ln (??) and P^'^ (r/) are the Laguerre and Jacobi polynomials respectively. 

The multi- indexed Laguerre and Jacobi orthogonal polynomials are given by (lA.llj) . 



2.2 Discrete quantum mechanics with pure imaginary shifts 

The Hamiltonian of the discrete quantum mechanics with pure imaginary shifts is (7 G 



U = y/V(x) e^^V*(x) + y/V*(x) e- lp ^V(x) -V(x)- V*(x) = A ] A, (2.19) 
A = i(e&y/V^-e-&y/V(a)), A ] = -i( ^/V{x)e^ - v^e"^) . (2.20) 

The ^-operation on an analytic function f(x) = ^2 n a n x n ( a n G C) is defined by f*(x) = 
^ n a*x n , in which a* is the complex conjugation of a n . The eigenfunctions <p n (x) and 
virtual state wavefunctions <fi v (x) can be chosen 'real', (f)* n {x) = <p n (x) and <p*(x) = 4> v (x). 
The intertwining operators Adi...d s an d A d da are given by 

M... ds = i (e^Vl„ ds (x) - e-fry/ V dl ... ds (x) ) , 

A... ds = -i(yJv*,.J.(x) & - ^Vl. As (x) e-i"), (2.21) 
V dl ...M = yJv(x-i^-y)V*{x-i^) 

x W 7 fc 1 ,...,0 da _ 1 ](g + ^) W 7 [0 dl ,...> d J(x-z7) ^ 
W 7 [0 dl ,..., ( /) da _ 1 ](x-z|) W 7 [0 dl , . . . ,(f) ds ](x) 

and the eigenfunctions (j)di...d s n{ x ) an d the seed solutions 4>di...d 3 \( x ) are expressed by 

4> dl ...d s n(x) = A(x)W^[(f) dl , . . . ,4> ds An]( x ), 

i>di-d s v(x) = A(x)W 7 [<j) dl , ■ ■ • 5 0d s ,0v](z), (2.23) 



A(x 



ii)w 7 [0 dl ,...,0 d j(x + zf; 



W 7 [0 dl , • • • , (f>d s ](x - «i)W 7 [0 dl , . . . , (j) ds ](x + i\ 
Here W 7 [/i, f2, ■ ■ ■ , fn]( x ) is the Casoratian 

W 7 [f u ...J n ](x)^U^^det(f k (xf)) , xf )d = f x + .(f -j) Tl (2.24) 
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and W 7 [-](x) = 1 for n = 0. We note that the deformed eigenf unctions and seed functions 

are 'real'; (j)* dl ... ds n( x ) = ^i.-.d^O) and 4> dl ...d s v( x ) = 4> dl ... da v(x). 

We consider two shape-invariant systems whose eigenfunctions are described by the Wil- 
son (W) and Askey- Wilson (AW) polynomials. Various data of these systems are: 

W: < £ < oo, 7 = 1, A = (ai, a 2 , a 3 , a 4 ), 5 = (5,5,5,5), k = 1, 

Tr/ ^ n i= i(ai + „ 7 . x , def 

V(x; A) = — — — , £ n (A) = n(n + b x - 1), 61 = a x + a 2 + a 3 + a 4 , 
2«x(2«a: + 1) 

0o(x;A) = y r(2 , x)r( _ 2 , x) > ¥>(*) = 2*. 

F n (x; A) = P n (j](x); A) = ^(^(x); a x , a 2 , a 3 , a 4 ) 

, w w , /-n, n + h-1, ai+ix, a!-ix 

= [at + a 2 ) n (ai + a 3 )„(ai + a 4 )„ 4 F 3 

V ai + a 2 , ai + a 3 , ai + a 4 



1 



27rn! (n + b x - l) n ]li<i<i<4 r ( n + a * + a i, 



ft »< A > ^ — r ( ^o (2 - 25) 

AW: < x < 7r, 7 = logg, g A = (ai, a 2 , a 3 , a 4 ), 5 = (5, 5,5,5), « = <T\ 
TT 4 fl — a e^) 

y(x;A)= /1 11J= ^ W1 3 ' £ n (\) = (q- n -l)(l-b 4 q n - 1 ), b, d = f a^a,, 



(1 - e 2ix ){l - qe 



&o(x;A) = J — 4 — - — L — , ?7(x) = cosx, ^(x) = 2sinx, 

F n (o;; A) = P n (r](x); A) = p n (r](x); oi, a 2 , a 3 , a 4 |g) 

g;g 



«i ™(oi«2, a 4 a 3 , aia 4 ; g)„ 4 ^ 3 1 

V fl 4 G£ 2 , (Z 4 (2 3 , Oxd 4 

ft (A) = 27r(6 4 g n " 1 ;g) n (6 4 g 2n ;g) 00 2g 

g)oo rii< 4 <j<4( a i a i9 n ; 9)00 ' 

where W n (r]; a x , a 2 , a 3 , a 4 ) and p„(^; ai, a 2 , a 3 , a 4 |g) are the Wilson and the Askey- Wilson 
polynomials jH] and q x stands for g( Al ' A2 ----) = (g Al ,g A2 , . . .) and < q < 1. The parameters 
are restricted by 

{a*, aj, ag, a 4 } = {a x , a 2 , o 3 , a 4 } (as a set); W: Rea^ > 0, AW: |o,| < 1. (2.27) 
The multi- indexed Wilson and Askey- Wilson orthogonal polynomials are given by (1A.19j) . 
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3 Recurrence Relations 



Ordinary orthogonal polynomials satisfy the three term recurrence relations [4"6] . 
7]P n (r]) = A n P n+1 (r]) + B n P n (r]) + C n P n _i(r?), 

or A n P n+1 (r)) + (B n -r ] )P n (r)) + C n P n „ 1 (r)) = 0, (3.1) 

with P-\{rf) == and A n C n+ i > (n > 0). When Po{f]) = constant is specified, the entire 
set of orthogonal polynomials is determined. We set 

A_a^0, P n ( V ) = (n<0). (3.2) 

Note that (13.11) holds for any integer n G Z, where A n (n < —2), S n (n < —1) and C n 
{n < 0) are arbitrary numbers, e.g. 0. We also set 

P dl ... dM Av) = (n<0), (j> dl ... dM n(x) = (n<0), (3.3) 

for the multi-indexed orthogonal polynomials. 

Corresponding to the three term recurrence relations, the multi-indexed orthogonal poly- 
nomials satisfy certain recurrence relations. We will show that the M-indexed orthogonal 
polynomials of Laguerre, Jacobi, Wilson and Askey- Wilson types satisfy 3 + 2M term recur- 
rence relations: 

Af+l 

E R l ^(v)P dl ... dM , n+ k(v)=0, (3.4) 

k=-M-l 

which holds for nGZ. Here the coefficients R^j(v) are polynomials of degree M + 1 — \k\ 
in 7] and their explicit forms are given by (13. 10[) for the Laguerre and Jacobi cases in § 13. 2[ 
and (13.171) and (I3.20j) for the Wilson and Askey- Wilson cases in § 13.31 These coefficients are 
expressed in terms of the coefficients of the three term recurrence relations, A n , B n and C n , 
as determined recursively by (13. 10[) or (J3TTTJ) . In other words, the coefficients R^kiv) are 
independent of the deformation data {0^ . . . ,0d A/ } except for M. As will be discussed in 
§3.4[ the data {(fid! ■ ■ ■ , (fi dM } are encoded into the initial data of the first M + 1 members of 
the M-indexed orthogonal polynomials. 

3.1 Three term recurrence relations 

Here we give the coefficients of the three term recurrence relations (13.11) for the standard 
Laguerre, Jacobi, Wilson and Askey- Wilson polynomials jH] with the input ^0(^7) = L 

L: A n = -(n + l), B n = 2n + g + \, C n = -{n + g-\), (3.5) 
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-r A 2(n+l)(n + g + h) _ ( h -g)(g + h-l) 



(2n + g + h)(2n + g + h + iy (2ra + g + h - l)(2n + # + h + 1) ' 
r 2(n + g-i)(n + ft-i) 

n (2n + <7 + /i-l)(2n+<? + /i)' 1 J 

w + 6i - 1 wlli^fc^ + aj + afc-l) 

n (2n + 6i-l)(2n + 6i)' n (2ra + h - 2)(2n + &i - 1) ' 

R = (n+h- i)nL 2 ( n + a i + ^) , n n 2 < J<fc <4( n + a i - jj 2 rs 7 n 

(2n + 6i-l)(2n + 6i) (2n + 6i - 2)(2n + h - 1) 0l ' 1 j 

AW : A n = — , „„_ 1W , z—z-r, C n - 



B„ 



2(1 - 6 4 g 2n - 1 )(l - M 2n ) ' 2(1 - fr 4 g 2n ~ 2 )(l - ft^ 2 "" 1 } 

2 2ai(l - 6 4 g 2n ~ 1 )(l - b A q 2n ) 



(3.8) 



2(l-6 4 g 2 "- 2 )(l-b 4 g 2 "- 1 ) 

3.2 Multi-indexed Laguerre and Jacobi polynomials 

In this subsection we derive the recurrence relations for the multi-indexed Laguerre and 

A 

dx 



Jacobi polynomials. First we note that the operator A = 4 — d x w{x) acts on a product of 



two functions f(x)<f>{x) as 

A(f(x)<P(x)) = f{x)M{x) + d x f{x) (f>(x). (3.9) 
Let us define Rn\(v) (n, k e Z, s <E 1>>-i) as follows: 
<to = (|fc|> a + l), <o 1 (^) = l, 

R [ : ] , k (v) = A n R [ :;l} k _M + (B n - v)R [ :; k 1] (v) + c n Rtl] k+ M is > 0). (3.10) 

For example, Rn\(v) f° r s = 0, 1 are 

S = 0: R [ °l( V ) = A n , nM (r,) = B n - V , F^_M = C n , 
s=l: ^(^ = 44+1, R [ ^ 1 (r ] ) = A n (B n + B n+1 -2r ] ), 
R%( v ) = A n C n+l + A n _ x C n + (B n - v ) 2 , 
RnMv) = C„C n _i, <Lx(?7) = C n (B n + 5 n _i - 277). 

This Rn\{v) (1^1 < s + 1) is a polynomial of degree s + 1 — \k\ in 77. By induction in s, we 
can show that 

d v R [ i k (v) = -(s + l)R [ Z%) («>0). (3.11) 
10 



We will show the 3 + 2s term recurrence relations of 4>n\x) == 4>di...d s n(x), 



E = s\(d x r,(x)) s ct> n (x) (s > 0), (3.12) 

k=—s 
s+1 

E <U^M+^) = ° (*>°)> ( 3 - 13 ) 

by induction in s (for n 6 Z). Since <$n{x) has the form (I2.12p . this (13 . 13[) means the 
recurrence relations of the multi- indexed orthogonal polynomials (|3.4|) . 
first step : For s = 0, (13.121) is trivial and (I3.13P is 

A n (f) n+1 (x) + (B n - r)(x))(f) n (x) + C„,0 n _i(a;) = 0, 

which is the three term recurrence relation itself. Therefore s = case holds. 

second step : Assume that (I3.12j) -( l3.13p hold till s (s > 0), we will show that they also hold 
for s + 1. 

By applying Aii...d s+ i to (13 . 1 3f) and using (13. lip , we obtain 

s+1 s+1 



0= E r[ :]Mx))&x)+ 9 Mx) E ^ R n]M x ))<t> l :U x ), 

k=—s—l k=— s — 1 

= E <U^))<^+^ 



fc=— s— 1 k=—s 

The second term can be expressed as 



E<* 1] W*M;*w 

fc=— s 

= i d ,.. ds ^ ^ i] (x) - ^(x) e 11 <fc 1] (*) 

k=—s k=—s 
s-1 

= S^77(X) ^ ^n,fe 21 ^n+k( X ) = s\{d x T}{x)) S (j) n {x) , 

k=-s+l 

where we have used the induction assumptions and (13.111) . Therefore we have 

s+1 

E R n]Mx))&x) = (* + l)!(d^(*)) S+ Vn(*), (3.14) 



fc=-s-l 
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which shows (13. 12j) with s — > s + 1. The three term recurrence relations and this (13.141) 
imply 

= (s + \)\{d x r](x)) s+l (A n <j) n+l (x) + (B n - r)(x))<fr n (x) + C„0 n _i(x)) 

fc=— s— 1 fc=— s — 1 

s+1 



k=-s-l 

s+2 

E (^ft^iW*)) + ( 5 » - ^))<U^)) + CnR^M*)))^ 

k=-s-2 



fc=-s-2 

which shows (13 . 1 3[) with s — > s + 1. This concludes the induction proof of (13. 12[) — (13. 131) . 

3.3 Mult i- indexed Wilson and Askey- Wilson polynomials 

In this subsection we derive the recurrence relations for the multi-indexed Wilson and Askey- 



Wilson polynomials. First we note that the operator A = i[e^ p \J V*(x) — e * p yV(x) ) acts 
on a product of two functions f(x)<p(x) as 

A(f(x)<f>(x)) 

= i(yjv*(x-q)f(x - q)<f>{x - q) - yjv{x + q)f{x + q)<t>{x + qfj 
= f w (x)i(Jv*(x-q) <j>{x - q) - y/v( x + q)(/>(x + qfj 
+ / H 0) (^v*{x-q)<j>{x - q) + ^/^\x~+q) 4>{x + q)) 
= f + \x)A(j){x) + f-\x) (yjv*(x-q)(/>(x - q) + ^v(x + q)<j)(x + qfj , (3.15) 

where f^'(x) are defined by 

f i+ \x) = \{f{x - q) + f( x + q)), f-\x) = ~(f(x - q) - f(x + q)). (3.16) 

Let us define R^\(x) (n, k <E Z, s G Z>-i) as follows: 

$£k(z) = AnR [ :;t-i( x + *D + ( R » - i( x - *§7))<; 1] (z + <?) 
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+ Cn&£l ] M1 (x + q) (*>0). (3.17) 

For example, R^ k (x) for s = 0, 1 are 

S = 0: R l °\(x) = A n , R l ^(x) = B n -r ] (x), R^x) = C n , 
s = l: Rl] 2 {x) = A n A n+1 , R^x) = A n (B n + B n+1 - rj(x - if) - rj(x + q)) , 
R [ n ] fi (x) = A n C n+1 + A n ^C n + (B n - rj(x - q)) (B n - r](x + q)) , 
<U(z) = C n C n - U <U(x) = C n (B n + B n _, - r)(x - q) - rj(x + q)) . 

By induction in s, we can show that 

R^t\x) = -lHx-i s ^l)-v{x + i s ^l))R^ ] {x) (s > 0). (3.18) 

By using this, we obtain the following expression of R^ k (x): 

Ri] k (x) = A n R^ 1 (x) + (B n - l( V (x - i§ 7 ) + r,{x + i§ 7 )))<* 11(+) (*) 

+ C n R [ ::l ] ^l(x) - \{ V {x - if 7 ) - V (x + qi)) 2 R [ Z 2 \x) (s > 0). (3.19) 

(For s = 0, the last term vanishes.) This recurrence relation with respect to s implies 
that R^ k (x) is a polynomial in r](x — if 7) + v( x + *f 7) an d {v( x — *f 7) ~ v( x + if 7)) 2 
(m — 0, 1, . . . , s). On the other hand the sinusoidal coordinate rj(x) satisfies 



T)(x - if 7) + 7}(x + if 7) 

rj(x - if 7)77(0; + if 7) 



2r](x) - \m 2 : W 
(g^r + q~^)r](x) : AW 

(^) + X) 2 :W 
rj^y+^qf - q -f)) 2 : AW ' 



which implies that any symmetric polynomial in rj(x — if 7) and r](x + if 7) is expressed as 
a polynomial in rj(x). Therefore we obtain 

R-n\( x ) = Rn\{v{x)) (\k\ < s + 1) : a polynomial of degree s + 1 — in i](x). (3.20) 

We remark that ^(x) = R^x). 

We will show the 3 + 2s term recurrence relations of <f>n\x) == ^...^n^), 

s+l 

E <l(^)^U(^) = (*>0)> (3-21) 

k=-s-l 
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by induction in s (for n 6 Z), Since </% (x) has the form (I2.12p . this means the recurrence 
relations of the multi-indexed orthogonal polynomials (13. 4p . 
first step : For s = 0, (EOT!) is 

A n (f> n+l (x) + (B n - T}(x))<t> n (x) + Cn^n-xix) = 0, 

which is the three term recurrence relation itself. Therefore s = case holds. 

second step : Assume that (13.211) holds till s (s > 0), we will show that it also holds for 
s + 1. Here we use simplified notation V^ s+l ^(x) = f Vd 1 ...d a+1 { x )- 
By applying A dl ... da+1 to (j3.2ip. we have 

o= E Cwfi 1 ^) 

fc=-«-l 

+ E $£t\x) (y/vw(x - q) ^ +h { x - q) + + q) + • 

fc=— s— 1 

By using ( 13. 18ft this is rewritten as 

s+1 

E Rn^W&x) = \H* ~ i§ fl) - V(* + ^7))CL s+1] (x), (3.22) 

fc=-s-l 

where 

fe=— s 

Then we have 

A»GK ] W + (X - \{ri{x - i*±* 7 ) + 77(x + z^ 7 )))Gt +1 ](x) + C n GH(x) 

= v^ [s+1] *(*-*i) [a. E ^wEU^ - *i) 

^ fc=-s 

+ - - + V(x + E ^ ~ *i) 

fc=— s 

+ a E - <§)) + (cc.) 

fc=-a ' 
I 8+1 / 

= y/vi^*(x-q) E UnR l :+lWx) 

k=-s-l ^ 
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k=-s-l ^ 

+ i fas - i^T) - ^ + ^7)))<; 1] (^)) + (cc.) 



s+l 



fe=-s-l 
- s+l 



+ Jy[^i](x + i2) E <l(* + *i)^+ fc (* + *I) 



fc=-s-l 

-(^(x - 1^7) - + ^7)) 

s 

X 

k=—s 



E ^(^(V^ 11 *^-^)^^ - *i) - (V^ [s+ % + *I)^+^ + ^)) 



= -5 - ^7) " V(x + ^7)) E ^(^K^)' (3-23) 

fc=— s 

where we have used induction assumption and (c.c.) represents complex conjugate. From 
(ET22]) and (ET23]) we obtain 

s 



fc=— s 

a* e ^s^^a*^) 

fc=-s~l 

-. s+l 

+ [B n - -(rj(x - *^ 7 ) + r,{x + ^7))) E ^ (+) 



k=-s-l 

s+l 



n 

k=-s-l 



namely, 

s+2 



/ 1 

0= E ^K ] (*)Un^ 
fc=-s-2 ^ 



s+2 



A;=-s-2 
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where we have used (13.191) . This shows (13.211) with s — > s + 1 . This concludes the induction 
proof of fl3T2TD . 



3.4 Initial data 

For ordinary orthogonal polynomials, the three term recurrence relations with an obvious ini- 
tial data Pq{t]) = constant (and P^iijf) = 0) determine the whole polynomials {P n (rj)} (n = 
0, 1, . . .). For the multi-indexed orthogonal polynomials {Pd 1 ...d M ,n(v)} (n = 0,1, . . .) with the 
derived 3+2M term recurrence relations, the initial data to determine the whole polynomials 
are the first M + 1 members of the polynomials, Pdi...d M ,o(??), Pd 1 ...d M ,i{v), Pd 1 ...d M ,M{rj). 
They are degree £,£ + 1, ...,£+ M polynomials in r\ and they are severely constrained by 
the input data of {0^ (x), . . . , 4>d M (x)}. 

The 3 + 2M term recurrence relations (13 Ah hold for any integer n G Z, which is classified 
into three cases (i) n < —M — 2, (ii) — M — 1 < n < — 1 and (iii) n > 0. For case (i), (13.41) 
is trivially satisfied because of Pd^.A^umij]) = (m <0). For case (ii), (13.41) is also trivially 
satisfied due to the fact 

R [ ^(rj) = (-M - 1 < n < -1, -n < k < M + 1), (3.24) 

which is a consequence of our choice A_i = 0. For case (iii), Pd 1 ...d M ,n+M+iiji) ls determined 
by (13. 4p and it is expressed by lower degree polynomials. As mentioned above, the polyno- 
mials Pdi...d M ,n{ r i) { n ^ M + 1) are determined by the 3 + 2M term recurrence relations (13.41) 
with M + 1 input data Pd 1 ...d M ,o{'n), Pd 1 ...d M ,i{v)- 1 ■ ■ ■ > Pd 1 ...d M ,M 

Note that the input data Pd 1 ...d M ,n(ji) (ft = 0, 1, ... , M) can also be calculated from the 
data of the lowest degree polynomial at each intermediate step, Pd 1 ...d a ,oi r i) (s = 0, 1, . . . , M). 
Since the 3 + 2M term recurrence relations (13. 4p are equivalent to 

M+1 

E R l ^{v(x))^d 1 ...d M n +k (x) = 0, (3.25) 

fc=-M-l 

giving the input data P dl ...d M ,n(v) ( n = 0,1,..., M) are equivalent to giving 4> dl ...d M n(x) 
(n = 0,1,..., M). If <p dl ...d a o(x) (s = 0,1,..., M) are given, (p dl ...d M n{x) (ft = 0, 1, ... , M) 
can be calculated in the following way. For s (s = 1,2, ...,M in tern), applying „ -( j s 
to (13.251) with (M, n) = (s — 1,0) gives 4>di...d 3 s{%) m terms of already known functions. 
Then applying A dl _ ds+1 to <j) dl ...d a s(x) gives 4> dl ...d s+1 s(x), and repeating this, and finally 
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we obtain ^...4^(4 Giving 0di...d s oO c ) is equivalent to giving Pd 1 ...d s ,o(v) ■ Note that 
Pdi...d a ,o{V'i ^) a ^di...d 3 (v'i ^ + ^)) which is a consequence of the shape-invariance [HJ ST]. 

4 Summary and Comments 

The multi-indexed orthogonal polynomials are a new kind of orthogonal polynomials sat- 
isfying second order differential or difference equations, whose degrees start from a certain 
positive integer I instead of 0, so that the constraints of Bochner's theorem are avoided. In 
this paper we have presented the recurrence relations of the multi-indexed orthogonal poly- 
nomials of Laguerre, Jacobi, Wilson and Askey- Wilson types. Corresponding to the three 
term recurrence relations of the ordinary orthogonal polynomials, the M-indexed orthogonal 
polynomials satisfy the 3 + 2M term recurrence relations (13. 4p . Their coefficients are ex- 
pressed in terms of those of the original three term recurrence relations. They are universal 
in the following sense; The derivation is based on (i) the three term recurrence relations (13. ip . 
(ii) the intertwining relations (I2.5p . (iii) the structure of the intertwining operators (I2.14p . 
f)2.2ip and (iv) the formats of the eigenf unctions (12.121) . The explicit expressions of the co- 
efficients of the three term recurrence relations (I3.5l) - (l3.8p and the explicit definitions of the 
multi-indexed orthogonal polynomials given in Appendix are not used. The multi-indexed 
orthogonal polynomials of Racah and g-Racah types [ID] are not discussed in this paper but 
the method is applicable to them, too. We leave this problem to interested readers. 

The 3 + 2M term recurrence relations need the initial data consisting of the first M + 1 
members of the polynomials. Namely, when the first M + 1 members of the polynomials are 
given as inputs, the other members of the polynomials are determined by the 3 + 2M term 
recurrence relations. For ordinary orthogonal polynomials (which start at degree 0), three 
term recurrence relations always hold and its converse is also true (Favard's theorem [16]); i.e. 
polynomials satisfying the three term recurrence relations become orthogonal polynomials 
(with respect to a certain inner product). It is an interesting challenge to formulate the 
converse of the 3 + 2M term recurrence relations. For example, in order that the polynomials 
determined by the 3+2M term recurrence relations become orthogonal polynomials or satisfy 
certain second order differential or difference equations, what conditions should be imposed 
on the first M + 1 members? 

The three term recurrence relations for the ordinary orthogonal polynomials are closely 
related to the closure relation between the Hamiltonian and the sinusoidal coordinate rj(x), 
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which leads to the canonical construction of the creation and annihilation operators a^, 
a<y+ '' } <Pn( x ) = A n (j) n+ i(x), a(~~>(f) n (x) = C n n _i(x) [U]. By transforming the original cre- 
ation/annihilation operators in terms of a series of intertwining operators Ad 1 ,„d a , -A^ d 3 ; 
the creation/annihilation operators of the systems of the multi-indexed orthogonal polyno- 
mials are obtained: 



(±) tfl ^1^2 _ _ 4-4 



a [M]( +)0 [M ](x) = An0 M (x)) a [M](-) [M] (x) = C^M (*). (4-1) 



It is interesting to see if the 3 + 2M term recurrence relations presented in this paper would 
lead to a generalized closure relation between the deformed Hamiltonian and the sinusoidal 
coordinate, and if it would give the above creation/annihilation operators. 

Recurrence relations for the exceptional (M = 1) Laguerre and Jacobi polynomials have 
been discussed in p2] in the context of bi-spectrality of orthogonal polynomials [19] . 
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A Definitions of the Multi-Indexed Orthogonal 
Polynomials of Laguerre, Jacobi, Wilson and 
Askey- Wilson Types 

For reader's convenience, we present the explicit definitions of the multi-indexed orthogonal 
polynomials of Laguerre and Jacobi types [8] and Wilson and Askey- Wilson types [H] . which 
are obtained by the method of virtual states deletion. 

There are two types of virtual states, type I and type II, which are derived by the 
discrete symmetries of the original Hamiltonian. We take the set of virtual states for deletion 
characterized by the degrees 

V = {d 1 ,..., d M } = {d\,..., d l Mv d l l dlj (M — Mi + Af n ), (A.l) 

and define 

A [Ml,Mnl = A + MiSi + M n <5 n . (A.2) 
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The eigenfunctions = (j>dx...d M n(%) — 4>Vn{x) of the deformed system = %di...d M 

= %x> have the following form: 

(j>vn(x) = ^ v (x)P Vin (rj(x)), (A.3) 

where Pv,n(ji) = Pd 1 ...d M ,n(v) is the multi-indexed orthogonal polynomial and the function 
^v(x) = ^di...d M { x ) * s expressed in terms of the ground state <fio(x) and the denominator 
polynomial ^vijf) = ^d 1 ...d M (ji)- The degrees of the denominator polynomial Hp (77) and the 
multi-indexed orthogonal polynomial Pv,n{v) are generically £ and £ + n, respectively, in 
which I is given by 

Mi M u 

£ = E 4 - 2 Ml(Ml " 1} + S d f ~ 2 MlI ( Mn - ^ + MlMn 

- 1 

= 22 d i - ~ M ( M - 1) + 2MiM„. (A.4) 
i=i 

A.l Multi-indexed Laguerre and Jacobi polynomials 

Two types of the virtual states are 

LI: #(x;A) = r^(iz;A), ^(77; A) = f P Y (-T]; A), 

5 ld = f l, £j(A) = -4(ff + v + J), (A.5) 
L2: ^(x;A) d = f v (x;t(A)), ^(77; A) ^ P Y ( m t(A)), 

t(A) d = f l-^, 5 Ild = f -l, £?(A) = -%-v-|), (A.6) 
Jl: fyx-Xj^^x-t'iX)), g(t ? ;A)^P v (i7;tf(A)) l 

tf(A)^0/,l-fc), 5 ld =(l,-l), £(A) = -4(s + v + f)(/»-v-§), (A.7) 
J2: ^A^^A)), £"(77; A) = f F v (?7; t n (A)) , 

t u (X) d = f (l-g,h), 5 Ild = f (-l,l), ^(A) = -407-v-J)(/» + v+|). (A.8) 

(We have changed the sign of <5i t n from those in [8].) The function ^v{x) in (1A.3I) is 

, Tr / \ Mil \\ 11 \\ def 0o(^; A Ml,Afn ) dcf / 2 : L , . v 

^(77(2;); A) [ -4 . J 

The denominator polynomial ^v{v) an d the multi-indexed orthogonal polynomial Pv,n(v) 
are defined by the following Wronskians: 

2© (77; A) = W[/ji, . . . , ii Mv vi, ■ ■ ■ , vm u ](v) 
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I / " — 

X 1 /i_„n(A/i+ 9 -I)A/„ / 1+ „WMi+/i-i)M, , ' (A. 10) 



-pD,n(^; A) = f W[/Xi, . . . , toft, V X , . . . , U Mu , P n ]{v) 



e -M I7? (Mi+g+\)M n . L 

X ^ (Afr+g+jpAfii / i+T? N (Mn+fe+DAf; _j ' (A.J I) 



A. 2 Mult i- indexed Wilson and Askey- Wilson polynomials 

Two types of the virtual states are 

type I: A) d = f v (x; t r (A)) , £{ m A) ^ P^; t\X)), g(s; A) ^ ^(x); A), 

t ! (A) = (1 — Ai, 1 — A 2 , A3, A 4 ), 5 1 = (— |, — |, |, |), 

_ / -(ai + a 2 -v- l)(a 3 + a 4 + v) :W . 
vl J ~ I -(1 - Oia 2 g- v - 1 )(l - a 3 a 4 ? v ) : AW ' 1 dj 

type II: 0"(x; A) = f v (x; t n (A)) , £"(77; A) = f P v (ry; t n (A)) , £"(x; A) = f £° (^(x); A) , 

t II (A) d =(A 1 ,A 2 ,l-A3,l-A 4 ), 5 l 

cii(\\ - / -(a3 + a4-v-l)(ai + a 2 + v) :W , , 

CvlAJ \ -(l-a 3 a 4 g- v - 1 )(l-aia 2 g v ) : AW ' K ' 

The function ^x>(x) in fjA.3[) is 

A) = (A.15) 



>2> 2' 2' 2/' 



Sx>(x - i\; A)Sd(x + z|; A) 
where a 1 (A) and a 11 (A) are 

Ql < A > = { Lr 1 : AW • a,I < A ) = { Lfl"' : AW • < A ' 16 > 

The denominator polynomial ^v{v) an d the multi-indexed orthogonal polynomial Pv, n (v) 
are defined by the following determinants: 



Sx>(x; A) = f Et>(t](x); A), Pv,n(x; A) = f Pv,n(v( x )'> A), (A. 17) 



j?(M) jS(Af) ^(Af) ^(Af) 

d 1 ' ' ' d l d 11 ' ' ' d 11 

"l "Mj "l u M n 
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2 ^ X )i 

Af-1 



: W 



' n^s^rij-Ji K--2-+^,a fc 
xn^njr 1 ^-^ 1 

]lfc=3,4 n*=i _1 arg^'(i+ 1 )(a feg -^ i e ia: , a^-^e"^; g).,- 
I x Ilfc=i 2 nSf 1 a^g^'+^g-^V*, a fe g-^V-; q) 3 : AW 



, (A.18) 



x ^|M(M+l) 



i?(M+l) 

X 4 



^(Af+l) ^(M+l) 
Y d 11 

"m t "i 



M 



IX I 



x n fe =i,2 nSK - f + ^, ^ 



B = < 



M 



— IX) 



2 ""13 
M ■ M 



jii 



: W 



rifc=3,4n^=i a fc J ?^ (j+1) ( a ^ A *e ix ,a k q >e fa ; 
I x Ilfc=i 2 IIS a- J g^^ +1 )(a fc g-f e^, a fe g"f e" 1 *; g), : AW 



(A.19) 



where 



(Xi M) ) J = rfixf^X, M)£(sf >; A), (1 < j < M), 

(n (M) ), = ^(4 M ^' M )^(4 M) ^) 5 



. ryu-y' ': X. M)r l j (xf a) ; A, M)P„( a ;f >; A), 



(A.20) 



and 



rjfcc!; ! • 



•; j ;A,M) = a I (A+(M-l)(5) - 



~K-|(M-1) I (M -l)2-(i-l)(M-i) 



x 



fc=l,2 



(A.21) 



: W 



rf(xi M) 



e HM + i-H) J](a fe g-^ i e-;g) J _ 1 (a fc g^e-;g) M - J : AW ' 

(A.22) 



fc=l,2 



M-l 



+ zx) J -_i(a fc - ^y 1 - ixjjif-j- 



: W 



x 



n (°* 

e i X (M + i-2j) "Q {akq -^ e ^ ^(^-¥6-; g) M _, : AW 



fe=3,4 



The auxiliary function i/?Af(x) is defined by 



M-2 



<p M (x) = (p(x) [A l ] ] [ (y>(x - ?|7)v?(x + %)) 



(A.23) 



fe=i 



and y^ot^) — V 9 !^) = 1 135] ■ Here [x] denotes the greatest integer not exceeding x. 
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